Hydrodynamic suppression of phase separation in active suspensions 



Suzanne M. Fielding 

Department of Physics, University of Durham, Science Laboratories, South Road, Durham. DH1 3LE 

(Dated: October 22, 2012) 

We simulate a suspension of active squirming disks over the full range of volume fractions from 
dilute to close packed, with full hydrodynamics in two spatial dimensions. By doing so we show that 
"motility induced phase separation" (MIPS), recently proposed to arise generically in active matter, 
is strongly suppressed by hydrodynamic interactions. We give an argument for why this should 
be the case, and support it with counterpart simulations of active Brownian disks in a parameter 
regime more closely akin to hydrodynamic suspensions than in previous studies. 



"Active matter" [l| comprises internal subunits that 
collectively drive the system far from Boltzmann equi- 
librium by each individually consuming energy. Biologi- 
cal examples include actively crosslinkcd polymeric fila- 
ments in the cell cytoskeleton cells grouped in living 
tissues Q ; suspensions of motile microorganisms 0, @ ; 
shoals of fish and flocks of birds Q. Non-biological ex- 
amples include vibrated granular monolayers [7H9(, and 
self-propelled synthetic colloidal particles floj [ . 

Distinct from the more familiar scenario of a passive 
complex fluid driven by (say) a global shear flow imposed 
at the system's boundaries, in active matter the driving 
out of equilibrium arises intrinsically in the active sub- 
units throughout the system's own bulk. Consequently 
active materials can spontaneously develop mcsoscopic or 
macroscopic mechanical stresses and deformations even 
without driving or loading from outside. Other ex- 
otic and generically emergent phenomena include swarm- 
ing, pattern formation, giant number fluctuations, non- 
equilibrium ordering and phase separation. (For a recent 
review see Ref. [l|.) These offer fascinating challenges to 
fluid dynamicist, rheologist and statistical physicist alike. 

Many active particles are elongated and so have an in- 
trinsic (steric) tendency to align with each other [IlTlloT ] . 
Activity mediated coupling between these orientational 
modes and fluctuations in the local number density then 
provides a generic mechanism for giant number fluctua- 
tions and phase separation 0, 0, EH , with the standard 
deviation AN in a subregion of material of mean number 
of particles N scaling as A^ a with a > 1/2. (For passive 
systems away from any transition a = 1/2.) 

Besides any such tendency for alignment, another 
generic mechanism for phase separation was recently put 
forward in the context of "run and tumble" particles, 
such as some species of motile bacteria 17|, ll8( . These 
swim in near straight-line runs at almost constant speed, 
between intermittent tumbles in which they suddenly 
randomize swim direction. The basic idea is that par- 
ticles (a) accumulate where they move more slowly and 
(b) move more slowly where crowded. Positive (a)-(b) 
feedback gives rise to "motility-induced phase separa- 
tion" (MIPS). This idea was recently extended analyt- 
ically to active Brownian particles [19(, consistent with 
simulations showing phase separation in active Brownian 



disks 20, 21 1 . These indeed lack any tendency for steric 
alignment; or for phase separation in equilibrium. 

To date these simulations lack hydrodynamic inter- 
actions, in which moving particles set up flow fields 
that influence their neighbours. But such interactions 
arise widely in active matter and this is fundamen- 
tally important because steady state properties in non- 
equilibrium systems depend strongly on dynamics (in 
contrast to equilibrium states, which depend only on the 
underlying free energy). This includes, crucially, the exis- 
tence (or otherwise) of activity-induced phase separation. 

The contribution of this work is to show that hydro- 
dynamic interactions in fact suppress MIPS, which ac- 
cordingly might not arise as generically in active matter 
as hitherto suggested, by simulating a suspension of ac- 
tive "squirming" disks [22|-|2J], with full hydrodynamics, 
over the full range of particle volume fractions from di- 
lute to close packed. (Being disks these lack any steric 
tendency to align, excluding a priori the first, non-MIPS 
mechanism for active phase separation discussed above. 
For simulations of active rods with hydrodynamics see 
Ref. [25[) We demonstrate carefully that hydrodynamics 
causes a key assumption of the (a) - (b) feedback mecha- 
nism of MIPS outlined above to fail. To support this, we 
further demonstrate suppression of MIPS in active Brow- 
nian disks in a parameter regime more closely akin to the 
squirmers, with MIPS recovered in the parameter regime 
explored previously for the Brownian particles [2(| . 
Models — For both squirming and Brownian dynamics 
we simulate P disks of radius R in a box of area L 2 . To 
prevent particle overlaps we impose a pairwisc repulsive 
force Fij = —f(a 3 — a 2 )dij, a = 2sR/ {dij — 2R) for inter- 
particle separation d^ < 2R(l + s). The effective volume 
fraction is then cj) = PnR 2 (\ + s) 2 /L 2 . 

The squirmers move through a Newtonian fluid of vis- 
cosity r\ by means of a prescribed tangential velocity 
Bi sin Oi + B2 sin 26^ round the disk edges. This gives a 
swim speed vq — B\j2 24j for a single disk undisturbed 
by any others, and an instantaneous swimming direction 
e-i = (cos#i,sin#i) for the ith disk. These swim direc- 
tions evolve over time due to hydrodynamic interactions 
between the disks, which are a priori unknown. Below 
we present data for j3 = B2/B1 = 0, but have verified 
suppression of phase separation for j3 = — > 00. 
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FIG. 1: Active Brownian particles. Top: average particle 
speed as a function of average particle volume fraction. Ex- 
pected spinodal shown by dotted line. Middle: scaled number 
fluctuations. Bottom: characteristic time r for reorientation 
of particle swim direction; and inter-particle collision time r c . 



The active Brownian disks [2(J 21, 26-28| have fi = 
voe.i + (J*Fij, again with swim speed vq for a single 
undisturbed particle. In contrast to the squirmers their 
swim directions are unaffected by interparticle interac- 
tions, having independently prescribed dynamics 9{ = 
rii(t), (m(t)nj(t')) = 2v r 5ij5(t — t'), regardless of the fre- 
quency or closeness of interparticle encounters. 

For both squirming and Brownian dynamics we simu- 
late P = 256 particles, heavily constrained by computa- 
tional cost for the squirmers (this P value being compara- 
ble with other squirmer studies [29| ). A much larger P is 
achievable for the Brownian disks, but we present results 
only for P = 256 to ensure as direct a comparison as pos- 
sible with the squirmers. We use units of length in which 
L = 1; time in which Do = 1; and mass in which 77 = 1 
(squirmers) or fi = 1 (Brownian). We set f = 1, s = 0.1 
to ensure the dimensionlcss parameter v / [if prescribing 
the small extent to which particles explore the interparti- 
cle repulsive potential is comparable to that in Ref . [2(J , 
giving almost hard disks. 

The remaining dimensionlcss parameters are then the 
volume fraction <fr = PirR 2 (l + s) 2 /L 2 ; and for the Brow- 
nian disks also the ratio 6 = Vq\P~/v t L of the prescribed 
time of decorrelation of swim direction and the charac- 
teristic time interval between particle collisions. 
Measured quantities — In what follows we report numer- 
ical results for the ensemble average particle swim speed 
v = -py J Q J^i \ v i(t)\i averaged over the large time in- 
terval t = — > T. (This speed differs slightly from 
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FIG. 2: Snapshots for active Brownian particles, correspond- 
ing to circles (a)-(d) in Fig. [1] 



the choice in Ref. [20j, extracted from the early-time 
ballistic regime of the mean squared particle displace- 
ment, which one can show analytically is equivalent to 
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We also report the scaled num- 
AN/N 1 ' 2 in a region of the sample 



Trio EiM*)l ! 
ber fluctuations 6 n 
with average number of particles N. 

We take the characteristic time t q for decorrelation of 
particle swim direction to be the time interval At = r Q for 
the function P{At) = J Q T dtj^i e 4 (t + At).e ; (i) to fall 
to 1/e. We find a reasonable measure of the time between 
inter-particle collisions to be r c = nR/(v — u goo d) i n 

which Wg 00 d = -ppp Jo ^2i v i(t)- e i{t) is a measure of the 
extent to which a particle manages actually to achieve 
its full velocity in its attempted swim direction. The 
degree to which this differs from the free swim speed 
is a measure of the degree to which scattering occurs. 
We have checked by direct observation that t c provides a 
reliable measure at low volume fraction. (At high volume 
fraction this comparison is more difficult to carry out, 
because particles slither round each other continuously.) 
Results — We explore first active Brownian disks in the 
regime where the particle swim direction is slow to decor- 
relate, C 3> 1- We do this to establish a point of contact 
with earlier simulation studies poj of active Brownian 
disks, which were performed in this regime, to provide a 
context in which to discuss our new results below. 

As seen in Fig. [1] (top), and correspondingly in Fig. 
3 of Ref. (though there with P = 10 4 ), the particle 
swim speed v decreases strongly with volume fraction 
<j>. Accordingly this system is a good candidate for MIPS 
via the (a)-(b) feedback mechanism discussed above. The 




N=64 



1 

0.8 
0.6 
0.4 
0.2 

0.6 

0.4 

0.2 




0.2 



0.4 



0.6 



0.8 



.oocfg^ 33 o% o ° o 
cP o 




(a)</> = 0.242 



(b)</> = 0.363 




0.5445 



(d)0 = 0.7865 



FIG. 3: Hydrodynamic squirmers. Top: average particle 
speed as a function of average particle volume fraction. Ex- 
pected spinodal shown by dotted line. Middle: scaled number 
fluctuations. Bottom: characteristic time r for reorientation 
of particle swim direction; and inter-particle collision time r c . 



expected spinodal, according to the instability criterion 
d\ogv/dlog(j) < — 1 of Ref. [13, EH, is located at the 
vertical dotted line in Fig. [1] 

Phase separation is indeed observed: see the top two 
snapshots of Fig. [2 and Fig. 1 of Ref. [2(|. This is re- 
flected also in enhanced number fluctuations: our plot of 
Sn (N) for each <j> (not shown) shows all the same features 
as Fig. 3 of Ref. [28[ , though obviously cuts off sooner at 
high TV due to our smaller number of simulated particles 
P = 256. Reporting in Fig. Q] (middle) the single value 
Sn=64, close to the peak of 5n(N) for this P, clearly 
indicates separation with binodal onset around (f> = 0.25. 

Consider now the hydrodynamic squirmers, our results 
for which are shown in Fig. [3] As for the Brownian disks 
the ensemble average swim speed declines strongly with 
particle volume fraction, suggesting that MIPS should 
again occur via the (a)-(b) feedback outlined above, with 
spinodal onset at the vertical dotted line according to the 



criterion d log v/d log <f> < — 1 of Ref. [17| . Remarkably, 
however, we find no evidence for bulk phase separation 
in the squirmers: see the snapshots of Fig. [1] and the 
correspondingly suppressed number fluctuations <5/v=64 in 
Fig. [3] (middle panel). (We do however find a tendency 
to form small string-like clusters at low volume fractions, 
as reported previously by other authors (29j.) 

Discussion - - To understand how hydrodynamics 
might cause this suppression of motility-induced phase 
separation (MIPS) we now revisit the original argument 
for MIPS, as first put forward in the context of run and 



FIG. 4: Snapshots for hydrodynamic squirmers, correspond- 
ing to circles (a)-(d) in Fig. [3] 



tumble dynamics and later generalised to other systems. 

Model run and tumble particles move in a series of 
straight line runs at a constant speed vq, when dilute, be- 
tween tumbles in which they rapidly reassign their swim 
direction. These tumbles occur randomly with a typical 
intertumble time interval a -1 . The argument for phase 
separation stemming from this motility is as follows, (a) 
It can be shown that the local particle volume fraction 
in any small region of fluid scales inversely with the local 
swim speed v. (Although intuitively reasonable, this is 
in fact a strongly non-equilibrium effect stemming from 
activity. For passive colloidal particles the volume frac- 
tion is uniform in the absence of an external field.) (b) 
It is assumed that between tumbles particles move more 
slowly in regions of high volume fraction, being impeded 
by crowding, rendering the swim-speed a decreasing func- 
tion of volume fraction v(<f>) < vq- Positive feedback 
between (a) and (b) gives phase separation, onset with 
spinodal instability from an initially homogeneous state 
when d log v/d log <p < — 1 . This argument extends to ac- 
tive Brownian particles of angular diffusivity v r , with an 
exact mapping a — > (d — \)v r in d dimensions |19| . 

Consider carefully part (b) of this argument. In or- 
der meaningfully to define an intertumble swim speed 
v((f>) < vo that is reduced by crowding, each particle 
must encounter many other particles between tumbles: 
the characteristic intertumble time a -1 must be large 
compared to the characteristic timescale r c between colli- 
sions. Transcribing this reasoning to Brownian particles, 
the angular diffusion time must be small compared to r c . 
Generalizing still further, we propose that regardless of 
the underlying microscopic dynamics - run and tumble, 
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Brownian or hydrodynamic - the decay time r Q of the 
swim-direction autocorrelation function must greatly ex- 
ceed the inter-particle collision timescale r c if part (b) of 
the argument is to hold and MIPS is to occur. 

With this reasoning in mind consider again Fig. [T] for 
Brownian disks in the regime f ^> 1 discussed above. The 
bottom panel of this figure shows the decay time r of the 
autocorrelation function of particle swim direction, to- 
gether with the typical time r c between particle collisions. 
Consistent with the externally imposed angular diffusiv- 
ity u T being small in these runs (large ( = 15), r D is a 
relatively large constant across all 4>: except in the very 
dilute limit each particle encounters many others during 
the time it takes to alter its swim direction: r c <C t q . 
Between angular reorientation events, then, each parti- 
cle properly samples the ensemble average reduced swim 
speed v{(j)) < vo: part (b) of the feedback loop holds 
and MIPS can indeed occur, as observed. Indeed, this 
mean-field nature of v((f>) was noted in Ref. 

In contrast, for the squirmers we are not at liberty to 
prescribe from the outset the timescale r Q of decorrelation 
of particle swim direction: instead this emerges naturally 
as a result of hydrodynamic interactions between the par- 
ticles, which are a priori unknown. (Put differently, there 
is no externally tunable £ for these particles.) Indeed, 
two-squirmer studies show that with hydrodynamics each 
scattering event ("collision") typically results in an O(l) 
change in swim-direction for each particle involved. Ac- 
cordingly, in these many-squirmer simulations we expect 
r w t c . This is indeed observed: Fig. [3] (bottom panel). 
Squirmers will therefore be unable properly to sample the 
reduced v((f>) < vq between reorientation events: part (b) 
of the feedback argument fails, and MIPS is suppressed. 

To support this argument we finally revisit the Brow- 
nian disks, but now imposing a smaller £ = 1 so that 
the particles reorient their swim directions much more 
quickly than in the large £ case discussed above, giving 
dynamics more closely akin to the squirmers (t d ss t c , 
Fig. [U bottom). Phase separation is then indeed strongly 
suppressed: see the bottom two snapshots of Fig. [21 and 
the significantly reduced number fluctuations of Fig. Q] 
(middle panel). This is despite the ensemble average ve- 
locity v((j>) still decreasing strongly enough for this £ = 1 
to satisfy the condition for spinodal onset. Put simply, 
particles that rapidly reorient do not have time to sample 
this mean field v{<p) in between reorientation events. 

Despite the absence of true bulk phase separation for 
the squirmers, and for the Brownian disks with r ~ r c , 
some particle clustering is nonetheless still apparent in 
Figs. [2k, d and Fig. [4ji-d. Whether this clustering can be 
interpreted in terms of a nearby but suppressed MIPS, 
which could therefore be a generic feature of active mat- 
ter with hydrodynamics, remains an open question. 
Conclusions — We have simulated a suspension of ac- 
tive squirming disks across the full range of volume frac- 
tions from dilute to close packed, with full hydrodynam- 



ics. In doing so we have shown that hydrodynamic inter- 
actions strongly suppress motility induced phase separa- 
tion. These findings should apply generically to active 
systems in which hydrodynamics are important, and in 
which the effective particle velocity depends on volume 
fraction via collisions impeding particle motion. Obvi- 
ous exceptions include systems in which v depends on 
</> instead by chemically mediated mechanisms (such as 
quorum sensing [1(3]), allowing MIPS to arise at relatively 
low volume fractions, and/or systems in which the par- 
ticles reside in (or on the surface of) a gel [3l|. Such ex- 
ceptions notwithstanding, the mechanism proposed here 
is expected to arise widely in active matter, and partic- 
ularly in active colloids, which form the focus of intense 
current experimental interest (lo| . 
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